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Abstract 

The aim of this paper is to analyze a class of random motions which models the motion of a particle 
on the real line with random velocity and subject to the action of the friction. The speed randomly 
changes when a Poissonian event occurs. We study the characteristic and the moment generating func- 
tion of the position reached by the particle at time t > 0. We are able to derive the explicit probability 
distributions in few cases for which discuss the connections with the random flights. The moments are 
also widely analyzed. 

For the random motions having an explicit density law, further interesting probabilistic interpreta- 
tions emerge if we deal with them varying up a random time. Essentially, we consider two different type 
of random times, namely Bessel and Gamma times, which contain, as particular cases, some important 
probability distributions (e.g. Gaussian, Exponential). In particular, for the random processes built 
by means of these compositions, we derive the probability distributions fixed the number of Poisson 
events. 

Some remarks on the possible extensions to the random motions in higher spaces are proposed. We 
focus our attention on the persistent planar random motion. 

Keywords: Bessel process, Gamma process, iterated Brownian motion, Laplace distribution, Struve 
function, random flight, random time, telegraph process. 



1 Introduction 

Diffusion processes have a central position in the theory of probabihty. Nevertheless, their main short- 
coming is the tinbotindeness of the first variation. For this reason a diffusion process is often not suitable 
to describe the real motion and many researchers have proposed alternative models having finite speed. 

The prototype of the random motions with finite velocity is the telegraph process. By assuming that 
the change of direction is governed by a homogeneous Poisson process N{t) with rate A > 0, we can 
define the telegraph process as follows 

iV(t) + l 

m = v{o) Yl {sj - sj-i){-iy-\ (1.1) 

where 1^(0) is the initial velocity assuming the values +c or — c with probability | and the times sj are the 
instants in which the j-th Poisson event occurs. Furthermore, T{t) is finked with the hyperbolic partial 
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differential equations, because its density law is the fundamental solution of the equation 



The telegraph process has been studied by several authors; for example Orsingher (1990), Foong and 
Kanno (1994), Di Crescenzo (2001), Stadje and Zacks (2004) and Zacks (2004). This model seems to be 
suitable to describe the real motion and it emerges in different fields. In the following, we provide a brief 
review of the possible applications of the telegraph process and its generalizations. 

• Physics. In the physical mathematics the connection between the telegraph process and the elec- 
tromagnetic theory strongly emerges. In particular, the equation ( |1.2[ ) describes the propagation of 
a damped wave along a wire. Weiss (2002) provided an interesting review of the physical applica- 
tions of the process T{t). 

• Biology. Models governed by hyperbolic differential equations and in particular the telegraph 
equation have been used to describe the movement of chemotaxis (see Hillen et al, 2000). 

• Ecology. The telegraph process is useful in ecology to model the displacement of wild animals on 
the soil (see Holmes et al., 1994). In fact, this model preserves the property of animals to move at 
finite velocity along the same direction. 

• Finance. Di Masi et al. (1994) proposed to model the volatility of financial markets in terms of 
T{t) . Di Crescenzo and Pellerey (2002) introduced the geometric telegraph process as a model to 
describe the dynamics of the price of risky assets, i.e. the authors replaced the standard Brownian 
motion with the standard telegraph process. Ratanov (2007, 2008) proposed to model financial 
markets using the telegraph process with two different velocities (as the risky asset tends upward 
or downward) and jumps occurring at switching velocities. 

• Actuarial Sciences. Mazza and Rulliere (2004) established a link between hitting times associated 
with the risk process (time of ruin of the insurance company) and the telegrapher's motion. 

A statistical analysis of the random model T{t) has been performed by De Gregorio and lacus (2008) 
and lacus and Yoshida (2009), when the sample path is observed at discrete times. 

In this paper we will analyze a one-dimensional random motion which in somehow generaUzed T{t). 
At time t > 0, the random speed of the motion is defined hy v = c cos 9, where c is a positive constant 
and ^? is a random variable with density law given by 

So, we consider a particle starting from the origin, choosing initially a velocity c cos 6i with probability 
law given by fiy{0). The particle travels maintaining its motion with the same velocity until a Poisson 
event happens. Now, the particle changes velocity independently to the previous one according to fiy{0) 
again and so on. At time t, we indicate the particle position on the real line with X^{t). 

The function fu{d) allows us to define a random motion in which the small displacements have bigger 
probability than large ones. Therefore, we have the physical phenomenon called inertia. Indeed, every 
object or particle moving on a surface suffers an effect due to the friction of the surface itself. Then, the 
particle will tend to go away from the starting point slowly. If z/ = 0, one has that fo{9) = - becomes 
the uniform distribution on the semicircle with radius one. In this last case the particle moves without 
inertia. 

In the same spirit of Orsingher and De Gregorio (2007b), we study the conditional characteristic 
and moment generating function of X^{t). We are able to derive the explicit probability distribution. 
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conditioned to number of Poisson event, in few cases: = and v = 1. Therefore, let N{t) denote the 
imderlying Poisson process governing the changes of the velocity, we have that 



p{Xo{t) e dx\Nit) = n) _ r (t) r (f + 1) / 2 



dx 27rr(n) \ct 



{cH'-x'Y 



P{X,{t)Gdx\N{t)=n) ^ r(n+l)r(n + 2) / 2 ^ _ . 

dx 27rr(2n + 2) \ct J ^ ' ' 

with n > 1. These results, permit us to put in light the relationship between Xo{t) and Xi{t) and 
the random flights studied by Orsingher and De Gregorio (2007b). A random flight is a continuous 
time random walk defined similarly to X^{t), but with direction chosen uniformly on an hypersphere. 
By means of the above probabilities, we can claim that, in distribution, Xo{t) and Xi{t) correspond 
respectively to the projection onto real axis of a planar and four dimensional random flight. 

In Section 3 we derive the first two moments of X,y{t), while for = 0, 1, we are able to explicit the 
moments of order p by means of special functions. Moreover, we will point out some connections with 
the related random motions on hyperbolic spaces. 

In the second part of this paper, we focus our attention on the random motions Xo{t) and Xi{t) 
evolving up a random time, leading to interesting interpretations of the related conditional density laws. 
In other words, we will introduce families of random times, containing some important random variables. 
In the probabilistic literature there are several papers devoted to analyze the properties of stochastic 
processes with random times. For example, the Brownian motion with Brownian time (iterated Brownian 
motion) has been studied by Burdzy (1993), Khoshnevisan and Lewis (1996), AUouba (2002), DeBlassie 
(2004) and Nane (2006). A link between the solution of fractional partial differential equations and the 
iterated Brownian motion has been extensively investigated by Orsingher and Beghin (2009). The iterated 
Brownian motion has been proposed as a model for a diffusion in a crack (see Burdzy and Khoshnevisan, 
1998). Beghin and Orsingher (2009) have studied a planar random motion with Brownian times; the 
authors have provided the conditional probability on the number of the events of a fractional Poisson 
process. 

It is interesting to consider random times derived by Brownian motion. For example, the Bessel 
process R'^{t) = where Bj are independent Brownian motions. Under the condition 

Vm> ^ — 1, m = 0, 1, the following result holds 

fill* I d P 2c-^ tU! 

P {x^iR'm e ^|«(.) ^ „} = ^(.^„^_,^.) ^^-i" 

where vq = ^, vi = n + 1, n > 1 and B{a, h) is a Beta function of parameters a and b. Hence, the 
random motion X^ stopped at Bessel random time changes drastically its probability distribution which 
becomes a Gaussian with variance given by a Beta random variable (up the scale factor (?t). These results 
can be generahzed by considering a Z-times interated Bessel process 7?.f (t) = i?f (i?2(---(^f+i(i)- -)))- 

For d= 1, the Bessel process becomes a reflected (around zero) Brownian motion |-B(i)|, that is the 
Brownian time arising in the iterated Brownian motion. Then, the above condition is always satisfied and 
the following distribtional equaUty emerges 

XU\Bm=B(^^Xl{t) 

For (i = 2, we obtain a time distributed as a Rayleigh random variable and Vm is strictly positive for each 
n > 1. This last fact permits us to provide the density law of 

Other relationships will be point out considering the composition with the sojourn time of a Brownian 
motion on the positive axis. 
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Let Ga {t) be a Gamma process with parameters a > 0, that is a random process with probability 
law given by a ga{s, t) = rX^J)^""^^"**' s > 0,t > 0. The Gamma random variable arising in various 
applications and it is useful to model the lifetime of a phenomena. For this reason, we deal with a Gamma 
time and for Xm{Ga{t)) obtain that 

P{Xr,,{G^{t)) edx\N{t)=n} 

dww^ (1 - w) " 2 ' _ 



r(5d^)B(f,v,„- f + 1) 7o y/¥lBc\2cy/wJ ' \Cy/^ 

with w,„ > ^ — 1. As for the Bessel process, this result can be extended considering the /-iterated Gamma 
random times gi{t) = Gai{Ga2{- ■ ■ (Gai+dt)) ■ ■ •))■ 

For a — 1, Gi{t) becomes an exponential process and X„i{Gi{t)) is distributed as a Laplace random 
variable with parameter jx^t)]' 

To complete the discussion on the random times, we deal with a clock obtained mixing ^^{t) and 
Ga{t) and studying the effect on the probabilistic structure of the random motion X,„(i). 

It is not an hard task to extend the previous results to the planar and four dimensional random flights 
with randomly varying time. Moreover, in the last Section, we will discuss the possibility to consider a 
random flight with drift, that is persistent along a portion of the surface. 



2 Moving randomly with friction 

Let us consider a random motion which describes the displacements of a particle starting from the origin 
of the real axis. The particle moves forward or backward with random velocity v = ccos 9, where c is a 
positive constant, while 6* is a random variable having density law 

Md)^ /^|r+^^\. sin^-^, (2.1) 

with 9 S (0,7r) and v > 0. Therefore, the particle moves with a velocity, randomly chosen on x- 
component of the unit semicircle according to ( |2.1| i, and it performs its motion until a Poisson event 
happens when another velocity will be chosen independently from the previous one. The position at time 
t of the particle is defined as follows 

JV(t) + l 

X^{t) = c ^ {sj — Sj-i) cos9j (2.2) 

where N{t) represents the underlying homogenous Poisson process with rate A > governing the 
changes of velocity, Sj,j — l,...,N{t) + l(so = 0,SN(^t)+i — t), is the time of occurrence of the 
j-th Poisson event and 9j's are independent random variables distributed as in ( |2.1| i. Furthermore, also 
N{t) and 9j are independent. From ( |2.2| l emerges that X^{t) is a telegraph-type process similar to T{t) 
defined in We note that X^, (t) has an infinite number of possible velocities and it has no necessarily 
alternating directions. Further, the particle at time t is located inside the interval (— ct, ct), and Xi,{t) 
has a fully absolutely continuous probability distribution, whilst in the law of T{t) a singular component 
appears (see Orsingher, 1990). 

We underline that for values of 9 close to | the density law (0) assigns probability mass greater than 
ones near or tt. This means that the process moves away slowly from the starting point. This represents 
the effect of the friction of the surface on which the particle performs its motion. Wlien v assumes high 
values, the density f,^{9) is highly concentrated aroung ^ (see Figure 1) and then the motion is slowed 
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Figure 1: The behavior of for v = 0,1, 2, 3. 

down. For this reason i' represents the level of friction to which is subject the motion. In other words, 
Xi, (t) defines a whole class of random motions indexed by the parameter namely the level of inertia. 
For ly — 0,we reobtain the uniform distribution on the semicircle with radius one and Xo(t) is exactly the 
x-component of the planar random flights studied in Orsingher and De Gregorio (2007b) or equivalently 
the projection onto real line of the sample path of a planar random flight. 

Our first result concerns the characteristic function of X^{t) conditioned on the number of Poisson 
events during the time interval [0,t]. 

Theorem 2.1. The conditional charade rsitic function of X,j{t) is equal to 

^{e--<'.|„,„ = „} = !il(2^r(.+ . . ...n (2.3) 

for n > I, while for n — 0, one has 

i?{e*"^'^(*)|iV(t) = O} = (^^^ T{iy+l)Mact), (2.4) 

where Ju{x) — X]fc^o(~'-^)^ T(k+()T(k+v+i) ^^^^ well-known Bessel function. 

Proof. In order to prove p.3| l and ( |2.4[ ), we observe that the Bessel function J^{x) admits the following 



J^z) = . — -r^ I e""-"""^ sin^" (j)d(l> (2.5) 



integral representation 
with Re{y) > 0. For n > l,we get that 



— / dsi--- ds,,—- — / sin " 6li<i6ii — — — / sin 6'„+id6l„+i 

{n+l 
iac ^^(sj — Sj_i) cos 6*^ 



i Jo Js„^^ fj[ {ac{Sj - Sj^i)y 



where in the last step we have used the integral representation ^2.5\ . For N(t) — 0, the position of the 
particle at time t is X{t) — ct cos 9 and then 



E 



□ 

It is interesting to observe that the charecterstic function of has the same structure of the 

one emerging in the problem of d-dimensional random flights (see formula (2.3) in Orsingher and De 
Gregorio, 2007b), where the parameter f is replaced by | — 1. 

For the moment generating function we present the following Theorem. 

Theorem 2.2. The conditional moment generating function of X^(t) becomes 

E {e^--<..|»„) . „} . ^^n. + 1))- /' .... - ... n a« 

given n> 1, while ifn — Q the following expression yields 

= o} = (^A^ r(i.+ l)/.(/3ci) (2.7) 

where Ii,{x) = X^fcLo r(fc+i')r(fc+t/+i) ^^P^^^^nts the modified Bessel function. 



Proof. The proof for ( |2.7| i and ( |2.6| l follows analogously to the one developed for ( |2.3| l and ( |2.4| l, noticing 
that 

2)^(2) Jo 

□ 
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The random motions obtained by setting u ~ and = 1 in the density law f,y{0) have a special role 
in this paper. Indeed, for Xo{t) and Xi{t), we are able to explicit in closed form their characteristic and 
moment generating functions and successively the density laws. In order to distinguish these important 
particular cases from the general random model X^{t), > 0, we will indicate them in the rest of paper 
with Xm{t),m = 0, 1. Moreover, we will use the following notation; vq — ^ and vi — n + 1. 

Corollary 2.1. For Xm(t), m = 0, 1, and n > I, we have that 

E[e^'^^-i')\N{t) ^ n] = ^ ^""'^^^yf^ J^^io^ct), (2.8) 

E[eP^^i^)\N{t) = n] = ^ ^"'^^'1^^^ I^APct). (2.9) 

Proof. We only give some sketches of the proof. Starting from ( |2.3| l, we prove ( |2.8[ ) for m = 0. It is 
possible to use the same approach used by Orsingher and De Gregorio (2007b), noticing that the n-fold 
integral 

ri+l 



/ dsi - ■■ dsn Y]_ Jo{ac{sj ~ Sj-i)) 

"'0 Js„-i . ^ 



can be worked out by applying recursively the following formula (see Gradshteyn and Ryzhik, 1980, 
formula 6.533(2)) 

x'^ia - xYJf,{x)J,{a - x)dx ^ ^(^+ ^^^^^ ^ J ,{a) (2.10) 

V27rr(/i + + 1) ^ 

with _Re(/i) > — Re{v) > — ^. Analogously, for m 1, by taking into account the formula (see 
Gradshteyn and Ryzhik, 1980, formula 6.581(3)) 

'^Mx)Ma^x)^^^(i^iy_^^ ^^^^^^ 



X a — X \^ V J a 

with Re{ii) > 0, Re{v) > 0, it is possible to compute the exact value of the following quantity 

rt rt n+1 



dsi - ■ dsnY\_ 

Jo Js^-i .1 



Ji(ac(sj - Sj_i)) 



Hence, the result ( |2.8[ ) for vi = n + 1 emerges. 

For the moment generating function, we need to prove that the following semigroup-type property 
holds ^ ^ 

x^{a - xYl,{x)I,{a - x)dx = ^^t^J^^^'' ^ W (2.12) 

V27rl (/i + + 1) ^ 

Indeed, since (see Gradshteyn and Ryzhik, 1980) 

e^-x^^lJx)dx = ^ n^'^+2) ) 







we have that 



e^^da / x^'{a-xYl^{x)h{a-x)dx = / e'^''x''I^{x)dx / e'^yy" I^{y)dy 
Jo 



,r(M+i)r(i.+ i) 



(/32 - ly+i^+i 

^+\) r 

2TrT{fi + iy+l) Jo 
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Furthermore, being (see Gradshteyn and Ryzhik, 1980) 

with the similar steps used to obtain ( |2.12[ ), it is possible to show that 

' 'dx= [- + - ]a-^I^+^{a) (2.13) 



X a — X \/i ly 



In conclusion, by means of ( |2.12| i and ( |2.13| l, and the same considerations done for the proof of ( |2.8| l, the 
proof (|2.9[) of immediately follows. □ 



Remark 2.1. We are able to give an integral representation of the unconditional characteristic function 
of Xo{t). Indeed, we have that 



E 



. (At)" r(f + 1) r 

y \ / , / e^"^* ' sin" 0d0 

Onr2 ( n+l \ I 
n=0 ^ ^ \~) 



^ (Ai)^ r (2 r g.act cos e Q^Q 

^ III .A^rrn+i^ /„ 



e"^* |jo(aci) + g^actcose (/^(Atsin6') + Lo(Aisin6l)) sin6ld6'| 



where L^(x) = X^fcLo r(fc+'^'^)r(fc+;^+^) ''^ ''^^ modified Struve function. 
For Xi it) we get that 



y. (Ay (n+i-i) + i r 



e 



;t^o ^ r(n + |) 



where Ea.p{x) = X]fc°=o r(afc+/j) the Mittag-Leffler functon. 

As stated before Xo(i) and Xi(i) represent two important particular cases of the class of random 
motions X^{t). This is due to the fact that by means of \2.%\ we derive their probability distributions. 
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Theorem 2.3. The following conditional density laws hold 
w/f/z |a;| < ct and m — 0,1. 

Proof. We only prove ( |2.15[ ) for m = 0, because the other results follow by means of similar steps. 
Instead of inverting ( |2.8[ ), we show that the characteristic function of the probability distribution p^ix, t) 
corresponds to ( |2.8| l for m = 0. Therefore, we have that 

E 



r(f)r(t + i) /2V 



cos(Q!a;)(c^i^ - x^) "2 dx 



7rr(n) \ct , 

7rr(n) \ct) ^ (2fc)!io ^ ' 



r(f)r(t + i)2 



fc=0 



^r(?^) ' {2k) 



r(t)r(t + i)r(^)2"-i - (act)^^ r(fc+i) 
^r(n) ^ (2fc)! r(fc+f + i) 

r(f)r(f + i)r(n)0F - .^ K^)^" r(2fc)V^2i-^^ 
7rr(n)r(f) ^ (2fc)! r(fc)r(fc+f + i) 



r(i.i)EM) 
■ =0 

-Ji±{act) 



act\ 



2k 



fc=o r(fc + i)r(fc + § + i) 

r(t + i)2t 



(act) • 



□ 



Remark 2.2. Theorem 2.3 permits us to point out the existing connection between Xq^I), Xi(t) and 
the random flights. Indeed, by setting v = in f^{0), we reobtain the uniform distribution on a unit 
semicircle and (as expected) Xq^I) represents the projection onto the x-axis of a planar random flight. 
Then 

corresponds to the distribution (4.4) obtained Orsingher and De Gregorio (2007b). For v — \, the 
probability law 



2n+2 



is the same of the one obtained by means of the projection of a random flight in onto the real line (see 
(4.1c) in Orsingher and De Gregorio, 2007b). In other words, the shadow on M of a four-dimensional 
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random flight is perceived by an observer located on the real line, as a slowed down motion. Therefore, 
in distribution, we have the following equality 

N{t) + 1 

Xi {t) — c {sj — Sj-i) sin di,j sin 62J sin (j) 



where 



is uniformly distributed on the four-dimensional hypersphere. 



Remark 2.3. Theorem 2.3 says us that p^{x, t) — t) if and only ifv — Vm- This means thatX^it) 

(absence of inertia) and Xi{t) is equivalent in distribution to a random model representing a particle 
slowly moving with the same speed until t. 



The unconditional density laws of X, (<) and Xi (t) are given by the following expressions 



Ae 



-At 



2c ^ V 2c 



fe-i 



r2(fc±i) 



Ae 



-At 



2c 



/O ( Vc2i2-X2j +Lo f 



-At 



(2.16) 



p\x-t) 



cyA^S(A^'''' ^'V r(fc+i) 



fc=0 



fc+ 1 



5 \cH 



cH 



(2.17) 



with < ct. 

We point out that the first term in ( |2.16| l, that is Iq {j^\J c?t'^ — a::^), is equal to the one presents in the 
absolutely continuous component of the law of a telegraph process (see Orsingher, 1990), that is 



,-At 



2c 



A/o f -^/cH'^~x■^ \ + ^/o (^^ ^/ cH-^ - x-^ 



while the derivative with respect to the time of the Bessel function is replaced by the modified Struve 
function Lo {^\J (P't'^ — x^) (up the constant A). 



3 On the moments and some relationships with random motions on 
hyperboUc spaces 

In this section we analyze the moments of the random motion Xy(t\ In particular, we are able to provide 
the first two moments of X^, (t) by applying the results contained in Section 3 of Stadje and Zacks (2004). 
Therefore, fixed n > 0, for the mean value one has 



E{X^{t)\N{t) =n} = ctE{cose} 
_ ctT{iy+l) 



%/^r(i.+ i)7o 



COS e sin^'' 9d0 



cos 6 sin^'^ 9dd — I sin 6* cos 
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where in the last step we have used the well-known integral 

1 p(-a+l-)p('&+l-) 

/ sin" 9 cos'' 0d9^ 7; A , Re(a) > -1, Re(b) > -1. (3.1) 

Jo 2 r(^ + l) ' ^ ' ' ^ ' 

Remark 3.1. The mean value can also be derived from Theorem \2.1\ (or equivalently from Theorem \2.2\ . 
It is clear that 



1 



{ac{sj - Sj^i)y a=o 2T(i/ + 1) 

while it is not hard to prove that 

d J^{ac{sj — Sj^i)) k / ac{sj — Sj^i)\^''^'^ ^ c{sj — sj^i) 

da {ac{sj - Sj-i)y ~ klT{k + u + 1) \ 2 J {ac{sj - Sj^i))" 

which calculated at a — is equal to 0. Then 



da y i 



= (3.2) 

a=0 



For the the second moment we have that 



E {xUt)\N{t) ^ 72} = —ch^E{cos^e} 



n + 2 ^T{i^+l)J„ 



, .r(i. + i)r(|) 



and after some calculations 



n + 2 y/nT{iy + 2) 
(i/ + l)(n + 2) 



('^ + + 2) 



n=0 



e 



^{E,,2{Xt)~E,,3{Xt)} 



(^+1) 

(^^^^ ' ^-^1 (3.3) 



{iy+ 1) Xt \ Xt 

where in the last step we have used the following relationships: £'1 2(2^) = -^1,3 = j^—^- 

As expected, if 1/ increases the action of the friction is stronger and the value of E {X^{t)} tends to 
decrease. Indeed, for growing values of v the particle maintains itself close the starting point, so that the 
probability distribution of X^{t) will be less sparse. 
Further, it is not difficult to show that 



ct 



[0 p is odd 



E{XP{t)\m = 0} = / xPpUx, t)dx = I r(.+i)r(4i) , . (3.4) 

J-ct \ >r(§+Ui) (^^) pis even. 

For = 0, 1, we present the following result. 
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Theorem 3.1. The p-th moment of XQ{t) and Xi{t) are respectively given by 



E{xm^e-'^[^l^ ^ {ctyv {j^(A^) +L^(A^)} + ''^'j^jj^^'p (3.5) 

E{Xl{t)} = fP±^\ (ctr {E,.+i{Xt) - |^i,|+2(At)} (3.6) 



for p even, whilst E{X^{t)} = E{Xf{t)} = ifp is odd. 
Proof. Let p be even, we get that 

-At o° / \ \ '=^1 1 ret 



k=0 ^ ^ " V 2 

Now, we splitting the above sum in order to carry out separately the even and the odd elements. 



E{xi{t)} = ^(rf^TM^ y, , 



2^ V 2 ;| ^^V2; r(A; + l)r(fc+2±i + l) 

A\-' 1 /A\-^^/At\''= 1 



2c) v^r(| + i) V2cy ^V2; r(fc+i)r(fc + | + i) 



2c ^ ' \ 2 )\ ^^\2) r(fc + i)r(fc + 2±i + i) 



+ 



2c) ^ni + i)^\2c) ^^[2) r(fc+|)r(fc + £±i + |) 

At , . . , . -E±i 



-(ci)^r(^) ^ ) (7,±i(Ai) + L^(Ai)) + 
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For Xi{t) one has that 

^ '^'^ cVX^t'.WtJ T{k+l)Jo ^ ' 



(rf)p{ii;i,£+i(At)-|i?i,£+2(At)} 



2 

If p is odd, immediately follows that E{Xl {t) } ^ E{Xl{t)} = Q. □ 

Remark 3.2. From p.6| l we immediately reobtain the expression p.3| l /or i/ = 1, i>y setting p = 2. For 
E{XQ(t)}, starting from p.7[) anj Mi/ng f/ie duplication formula, one has that 



E{Xm - e~^\ctYY. 



^^^fc!(fc + 2) 



2c (^*)'^(2j^^V2c; r('^)r(| + 2) 

{ctY^ [xtf r(| + i) 
2 ^or(fc + i)r(| + 2) 



fc=0 



r(fc + i)(fc + 2) 



2V (A^)^(fc + 2-l) 
^ ^ ^^r(fc+l)(fc + l)(A: + 2) 

= e-^*(rf)2{^i^2(Ai)-^i,3(Ai)} 
corresponding to \3.3\ for v = Q. 

Remark 3.3. Since, for x > Q and ^ > 0, I^ix), L^(a;) fln^^ are monotone increasing functions as 
X oo, the following approximations hold as a; — > 



x^ 
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Therefore, for Xt — s- 0, we have that 



Xt , r(2±i; 



2 Ar(i + i) 



E{xm ~ (P±^) {cty- ^ 



V 2 ; ' ' r(| + 2) 

We conclude this Section discussing some connections between the previous results and the random 
motions moving on a Non-Euclidean plane. The upper half-plane H2 — {{x^y) : ?/ > 0,a; G M} 
endowed with the metric 

y 

is a model of Non-Euclidean (hyperbolic) space. The geodesic curves in this space are eiher the vertical 
half lines or half-circles whose centers lie on the x-axis. Similarly to Orsingher and De Gregorio (2007a), 
we consider a motion 

Y,{t) = e^-(*) (3.8) 

developing on the y-axis of the space iJ^, starting from the origin O — (0, 1) at time t ~ 0. The 
probability distribution of Ym{t) with m = 0, 1, is equal to 

p"(log2/,t)-, y>0 

y 



and by means of Corollary 2. 1 the conditional mean values of Yoit) and Yi (t) become 

E 



E 



{Yo{t)\N{t) = n} = £;{e^^°(*)|7V(t) ^ n} = £ii±yil/ (ct) 

L J (3=1 (ctj 2 ^ 

{Ymmt) = n} ^E{e^^^('^mt) ^ n} \^^^ ^ ?!^^^>±^/„+,(ct) 



Furthermore, by considering at time t the hyperbolic distance ri^{t) from the origin O of Y^{t), we have 
that (see Orsingher and De Gregorio, 2007a) 

.max(l,y„(t)) J,, 

vAt) = - = \XAt)\ (3.9) 

Jmin(l,y„(i)) y 

We are able to obtain a lower bound for the distribution function of rju{t). Indeed, we get that 



PM)<V) = P{\X,{t)\ < rj) 
1 

V 



^ l--E{X',it)} 

where E{X^{t)} is given by the formula ( |3.3| l. Clearly, for = 0, 1, we obtain the exact expression of 
the distribution function of the hyperbolic distance, namely 

Jo 

with < 77 < ci and m = 0, 1. 
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4 Randomly varying time stochastic motions 



So far, we have analyzed a random motion (t) evolving up to no-random time t > 0, deriving its exact 
probability distribution in two particular cases. In this Section, we focus our attention on the random 
motion Xm{t), m = 0, 1, defined as in ( |2.2[ ), with randomly varying time. 

In order to develop our analysis , we take into account families of random times which include some 
well-known random variables. In particular, we consider Bessel and Gamma processes as random times. 
Our choice falls on these two processes because at this way, we are able to include a wide range of 
probability distributions often used to model the time in many theoretical and real situations. Clearly, 
every random variable successively used as random clock will be supposed independent from X,n {t) and 
N{t). Therefore, we analyze the effect due to the composition of these random times with the random 
motion Xm{t) on the related density laws. 



4.1 Random times involving Brownian motions 

Let us consider a Bessel process starting from zero 



. t>0,d>l, 



where Bi{t)s are independent standard Brownian motions. It is well-known that the probability density 
law of R'^{t) is equal to 

1 r^"i r2 

/''(r)^— -— je-^, r>0. (4.1) 

r(|)23-it5 

At time t, we deal with a random motion X„i{t), m — 0,1, with a Bessel random time R'^{t). Since 
Xm{R'^{t)) is located inside {~cR'^{t), cR'^{t)), its support is the whole real line. Recalling that — ^, 

vi — n + 1 and by indicating with B{a, b) — '"p°j'j^ff , a > 0, 6 > a Gamma function and with B{t) a 
standard Brownian motion at time t, we are able to provide the following theorem. 

Theorem 4.1. Given N{t) — n, with n > 1, such that Vm > | — 1, we have the following conditional 
distribution 

1 

p{Xr.{R'^{t))edx\N{t)^n}^—^ ^— / u^-.-\l-wr--^^=^dw (4.2) 

with X € M. 

Proof. By using a similar approach to that adopted by Beghin and Orsingher (2009), we can write that 

P {X^iR'^it)) e dx\N{t) = n} 

poo 

= / P {Xm{R''{t)) e dx\N{t) =n,R'^{t) ^ s} P{R'^{t) eds} 



Therefore, by means of Corollary 2. 1 and for any n > 1 such that > | — 1, we are able to explicit the 
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Fourier transform of X,n{R'^{t)) as follows 



\acj r(f)22-42 Jo 



T{v.^ + 1) ^ (-1)^ / aVA ^ 

^ fc!r(fc + + 1) ^ 2 

B(i,,«™-i + i),t^ k\ y 2 J J, 



B{j,v^-'l + l) Jo 

Finally, by inverting E |e^"-^"(^''(*)) | A^(t) = n| the result ( |4!2| ) emerges. □ 

The probability ( |4.2| l claims that, conditionally on the number of Poisson events such that v„i > | — 1, 
the random process Xjn{R'^{t)), is distributed as a centered Gaussian with variance c^tW, where W ~ 



Remark 4.1. By using the same approach of the previous proof and bearing in mind the Theorem 2.1 
under the condition ly > ^ — 1, we get that 

P{X,{R'{t))edx\Nit) = 0} = --j-^ / w^-\l^wr-^-^^^dw (4.3) 



^ 2 ' 2 

In particular, we are interested to the random motions obtained by setting v = Q and v = I. It is clear 

d 
2 



that for i> = 0, the condition v > | — 1 /i satisfied only for d — \, whilst for d = 2 one has 



P{Xo{R^{t)) e dx\N{t) = 0} = dx^^. (4.4) 

\'2'Ktc 

If V — \, the above condition and the representation (|4.3[) hold for both d = \ and d = 2. 



From ( |4.1| i, we can derive some well-known probability distributions. Indeed, for d = 1, we get 
R^{t) = \B{t)\, that is a reflected Brownian motion around the x-axis and the its density law becomes 

/i(r) = ^=e~TF. Furthermore, \B{t)\ represents the Brownian time used in the definition of the 
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iterated Brownian motion (see, AUouba, 2002). For d = 2, we obtain a Rayleigh random variable with 

2 

/^(r) = |e~TF, which also emerges analyzing the distribution of the maximum of a Brownian bridge. 
Moreover, the probability ( |4.2| l holds for each n > 1, being the conditions i;,,, > — 5 (d = 1) and 
Vm > 0{d — 2) always satisfied. Actually, for the process Xi, the representation ( |4.2| l also yields when 
d — 3,4,5. Nevertheless, we restrict us to the cases d — 1,2. Therefore, we are able to explicit the 
unconditional probability distributions for X,n{\B{t)\) and Xm{R'^{t))- 

Theorem 4.2. For d = 1, the following probability yields 

Ct x^t 

P{X^{\Bm e dx} = dx r ^^t^ P"'{y,t)dy, xeR, (4.5) 

where p'^(y, f) and (y, t) are defined respectively by ^2.16) and ( |2.17[ ). Furthermore, for d — 2, one 
has 

1 

P {Xi{R^{t)) e dx) ^dx I e-^^'"[l + Xt{l - w)] dw (4.7) 

Jo V27riwc 

Proof. By taking into account ( |4.2[ ) and ( |4.3[ ), if m = (uq = ^) and d = 1, we get that 



^ Jo ^ ' V27Ttwc^^ n\ ^ > T{^) 



Ate-^* /"i _ie-^ ^ 1 /At \ 

/ w 2 > — - — — — — vl — w 

2 Vo V2TTtwc^^V^{'^)\2^ ) 



y„ 72^ „eor^(^)U^'^' 'J 

-ct ^2^ 2c i^,T^{^) ^2c^ 
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For m — 1 {vi — n + 1) and d 1, we have that 

1 ( \^^^ 1 — ^ — 
—P{Xi{\B{t)\)edx} = e-^'Y^ ^ w-'Hl-wT+'^^^^^dw 

e-'* r -le-^^iXtr^^ ^„+ir(n + 2)_^ 

= / w ^ , > ^ — —(l-wV^^—^ ^dw 

7o %/2^c^^ n! r(n+f) 

— At /-l ^>r— °° I 1 

_ / "5*^ n+1 (i 



Therefore, the result ( |4.5| l is proved. Developing the quantity 

(At) 



g-At ^ e dx\N{t) = n} 



n 

n=l 



for TO = 0, 1, and by taking into account the Remark 4. 1 it is not hard to prove the results (|4.6|l, (|4.7|i. □ 



We point out that the probability distribution ( |4.5[ ) says us that the process Xm{t) stopped at reflected 
Brownian time is equivalent in distribution to a Brownian motion with variance jX^^{t), i.e. 

x^mm = B (^\xi{t)^ (4.8) 

Remark 4.2. Recalling that for a centered Gaussian with variance the moments are given by 



r(p + i) 
2lr(§ + i) 



- 2 

if p is even, while are ifp is odd, we obtain that 



where E{XP^(t)} is defined as in Theorem 3.1 



The result ( |4.8[ ) is more general. Indeed, let M{t) be a Gaussian process with mean and variance 
cr^(t), it is no t hard to prove by using the same argument adopted in the proof of Theorem 4.1 and 

Corollary 4.2 that Xm(|A/'(t)|) is distributed as a Gaussian random variable with variance ^^-j^X^{t). 



In other words, the following distributional relationship holds 



xuwim^sf^xiit) 



For example, if: 
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Af{t) — (t), that is a fractional Brownian motion with Hurst index H e (0, 1), we have that 
cr2(t) = ^jjj jjjgjj x„(|B^(t)|) = B {t'^"-'^X^{t)), which contains as particular case the 
result ^ for H = i; 

A/'(<) = /o h{s)dB{s), where /i(s) is a well-defined deterministic function, the variance is cr^ (t) = 
/o /i2(s)ds and then /i(s)dB(s)|) = B (^^^^^^^^(t)) . Clearly, for h{s) = 1, we 

reobtain the equality ( |4.8| ); 

AA(<) = /g\B(s)ds, the variance is given by CT^ (i) = y, therefore we obtain that B(s)ds) — 



Baxlit)). 



Theorem |4T| can also be generalized by dealing with an ^-times interated Bessel process, namely 

7^f(^) = i?^(i^^(...«l(^)...))), t>o, 

where i?^s, j ~ 1, 2, Z, are independent Bessel processes. The random process Tlf{t) has density law 
given by 

(r(f)22 + 1 Jo Jo q q t2 

and leads to the next result. 

Theorem 4.3. Given N{t) — n, with n > 1, such that Vm > ^ — 1, we have that 

P {X„,(7^f (t)) e dx\N{t) = n] (4.9) 



i3(f,«™-^ + i)(r(|)2(i-i))'7o 4 7o tl Jo 



dti / ^ 5 rfi2 • • • / 5 dt 



1 . , 



X (1 — W) " 2 -rf^; 

Jo v27rtiwc 
w/f/i a; e M. 

Proof. Since Xm{T^f{t)) has support on the interval {—c}Zf{t),c}Zf{t)) and then on M, we can write 
that 

P{X„(7^f(^))e^^a;|iV(^)=.n} 

= p {x„^{nf{t)) e dx\N{t) = n, nf{t) = s} P{nf{t) e ds} 

Jo 

r(^;„ + i)r(w„) /""^ / 2 Y"'" 2 2 2Nt.„-i, d-ltd, 



277 r(2w^) 



2 \ / 2 2 2\t)„-i-, d-lrdf \i 

I — I (c s -x)"* n{|^|<cs}S fi[s)ds. 
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The conditional characteristic function becomes 



T{vm + l) I ( — ) J,,So^cs)ff{s)ds 



acs 



Tjv^ + 1) r r I TT ^^'-1 



J fc!r(fc + + 1) VYJ 











Jo 




POO 




/ dti-- 




Jo 





fe=0 



i=2 t^ I S "'-^ y"-^ ~ 2 



-e 



fc!r(«„, - 1 + 1) 



^y^rn + 1) / / J. I TT 



TK, - f + i)(r(f ))'+i2(i-i)' Jo Jo I t 



e 



with t. By inverting the so-obtained characteristic function the proof is completed. □ 

Similarly to the simple Bessel time, conditionally on N{t) = n, such that the constraint «,„ > 
^ — 1 is satisfied, XmiTZfit) is distributed as a Gaussian random variable with variance given by 

c^RfiRli- - - {Rf{t)) - - - ))W, W ^ B{j,vm - i + 1). 

Obviously, also the result contained in the Remark [4T| can be generalized as well. Indeed, by means 
of the same approach used in the proof of the Theorem 



4.3 if > I - 1, we obtain that 



P {X^(7^f (t)) e dx\N{t) = 0} (4.10) 

dx 1 t<}-^e-^ f^ti-^e-^ f^tf-^e-i 

' -dti / ^ 2 dt2 - -- ^ 2 dti 



B{l,^-i + i)iT{i)2(i-^)yJo tl Jo tl Jo ti 

[l-w) 2 -dw. 

^2'KtiWC 

For = and d = 2, it is easy to show that 

P{Xo(7^f(^)) e dx|7V(t) = 0} (4.11) 

, f°° tie'-^ , f°°t2e'^ , f°° tie-* , e ^ 
= dx dti / dt2 ■ ■ ■ dti 



t2 Jo ts Jq t v^27rt 
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Bysettingd = 1, 7?, becomes an /-iterated Brownian motion , namely 7?.;^ (i) = |i3i(|i?2(| • • • (|-B/+i(i)|) 
where Bjs are independent Brownian motions. Then, for each n > 1, we get that 

P {X„,{nl{t)) e dx\N{t) = 7i} (4.12) 

rdti / — at 2 ■ ■ ■ / — i^dti 



B{\,Vm + \) Jo Jq Jo \/7rt 



Furthermore, after some calculations similar to those of the proof of Theorem 4.2 we are able to explicit 
the following unconditonal distribution 

p{x,„(|Bi(|i?2(i • • • m+iim • • • )i)i)i) e dx} (4.13) 

dx2i I -^dh -^dt2--- -^dti ^ P"'{v,t)dv 

'7rr2 Jo v'''r3 Jo Y'l't J-ct ^/^-KtYy^ 



which allows us to claim that 

x„(7^K^)) = B (\B^(\ ■ ■ ■ (\Bm) ■ ■ • )l)l^) ■ 

Now, we analyze the effect due to the random clock defined as the time spent on the positive axis 
(sojourn time) by a standard Brownian motion B(t\ t > 0, namely 

T{t) = / l[B(s)>o}{s)ds. 
Jo 

The density function of r(t) is 7(5) = — . , < s < t, that it is also known as arcsin law. We have 
the following result concerning X,„(r(i)). 

Theorem 4.4. For n > 1, the following probability distribution holds 

dx 



P{Xm{T{t)) e dx\N{t) - n} 



Jo J sit - s) Jo 27rVc2s2^ 



B{1,v^)t: Jq ^s{t - s) 

Proof. In this case 

P{XUm) e dx\N{t) = n} = ^"^^.^^^5^ ^* ^ ~ xr-h{s)ds 
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As done so far, we consider the Fourier transform for Xm{T{t)). Hence, we get 

rK. + i)g (-1)'= '--2fe 

k=0 



k=0 

r(t>^ + l) y^ (-l)^r(2fc+i) /actV' 
^ ^ fc!r(A: + u„ + l)r(2fc+ 1) V 2 / 

(-i)'=r(2fc + i)r(fc + i)rK„)r(i) /actv 



TT ^ A:!r(fc+l)r(fc + w,„ + l)r(2/i: + l) V 2 / 

fc=0 

TT 



TT — (/c!)2 V 2' 2/ ' ' V 2 y 



Therefore, by taking into account the Theorem 2.ip.3 we obtain that 



dx vr Jo Jo 277 yJcH'^z'^w 

= (s = tz) 



(1 - 



Jq \/ s{t — s) Jo 27r ^/c^s'^w — a;2 

□ 

From Theorem ( |4.14[ ) we conclude that the random motions X„i{t) with random time T{t) is equiv- 
alent in distribution, to the random motion Xa{t) with a random time given by T{t)Vw, where W ~ 
B{l,Vm), which maintains the velocity initially chosen until t. Then, one has that 



P{XmiT{t)) e dx\N(t) = n} = P{Xo{T{t)VW) G dx\N{t) = 0}. 
Remark 4.3. The random process T{t) is also connected with Xra{R'^{t)). Indeed, recalling that 

P{T{t) e ds\B{t) > 0} = — ds, P{T{t) e ds\B{t) < 0} = —^^^^ds, 

T^t y/[t - S) TTf 



< s < t, from ( |4.2| l and ( |4.3[ ) the following equalities hold: 

P{Xo{\Bi{t)\) edx\N{t) = 0} = P{Bi{c^T{t)) e dx} 



dx 1 e 2c 



(is, 



Ti" Jo \/s{t — s) \/2ttsc 

P{Xo{R^(t)) e rfx|Ar(<) = 2} = e da;|7V(t) = 0} 

= P{Bi{c'^T{t)) e da:|B(t) > 0} 

dx2 /■* ^/s e"5^ , 
-as, 



Jo - s) V27rsc 
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P{Xo{\B,{t)\)edx\N{t) = 2} = P{Xi{\B,{t)\) &dx\N{t)=0} 

= P{Bi{c'^r{t)) G dx\B{t) < 0} 

dx2 /■* Vi- s e~^s 

= / — y= . as, 

T^t Jo Vs V27rsc 

where Bi{t) is a independent Brownian motion with respect to B{t). 



4.2 Compositions with Gamma random times 

In this part of the paper, we deal with a second class of random times different with respect to the previous 
one. We indicate with Ga (t) a Gamma random process, with parameter a > 0, governed by the density 
law ga{s, t) = rfs) > 0- Analogously to the Bessel case, first of all we study the conditional 

probability distribution. 

Theorem 4.5. Given N{t) = n > 1, such that > f — 1, the random process Xm{Ga{t)) has the 
following conditional probabilities 



P{Xm{G^{t)) e dx\N{t) = n} (4.14) 

dx , .V t ( M \^ J. ( t\x\ 

r(^^±i)i?(f,t;„- 1 + 1)70 V^c\2cVwJ 

where Kn{x) is the second type modified Bessel Junction. 

Proof. The random process Xm{Ga (i)) admits as support the real line, hence we get that 

P{Xm{Gc,{t))edx\N{t)=n} 

P{X„(G„(t)) e dx\N{t) = n, Gc.it) = s}P{G„(t) e ds} 

dx r(t.„ + i)r(i;„) ( 2 



27r T{2vm) r(a) Jq \cs 
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Then 



ac 



T{a) Jo s- 



+ 1) ^ 2^ _^ + 1) V 2 J 



r(a) ^ fc!r(fc + + 1) 

TjVm + 1) ^ (-l)'=r(2fc + a) /QC\ 2fc 

r(a) ^fc!r(fc + «,„ + i) V2J 

(by duplication formula) 

r(^;„, + 1) ^ (~i)^r(fc + f)r(fc + g±i)2^'--+°-i ^00^" 

V^r(«)r(..-t + i)g^(TJ n' + ~J''V^'+2'^'"-2+i 
r(t)r(^)rK„-f + i)^£^^^lTJ ^ 1,^ + — j ^ 2'"" " 2 + ^ 



1 



Jo Jo 



r(ii±i)i3(f,«,„-f + 1) Jo 
By inverting the above quantity, we obtain that 

j^P {Xm{G^{t)) e dx\N{t) = n} (4.15) 

= — — r- / dze z 2 — — / dww^ (1 - w) " ^ 

r(2±i)yo S(f,«m- f + 1) 7o Vi^c 



' dww^ (1 - w) " 2 — S= if-o ' 



a; 



r(^^) B(f - f + 1) Jo V7i^cV2c/wy " ^ VcVw 



2 ; 2 ' "™ 2 

where in the last step we have used the following integral representation 



K,i-)-U^Y I e-^-^^z-^-'dz. 



' 2 V2. 

Moreover, from ( |4.15[ ) it is easy to verify that 

P{X„(Ga(t)) € dx\N{t) = n} = 1. 

□ 

Remark 4.4. From \A.\5) it is straightforward to observe that 

P {X^{G^{t)) e dx\N{t) ^ n} 
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is equal to the density of a Gaussian random variable with variance -^ZW where Z is gamma variable 
with parameter and 1, while W ^ ,Vm ^ ^ + 1). 

Remark 4.5. analogously to the Bessel case, the approach developed in the previous proof is enought to 
show that 

P{X^{Ga{t)) e dx\N{t) = 0} (4.16) 

dx r\ ,,„^ t ( t\x\ Y f t\x 



r(^^±i)B(f,i^- f + 1) io ^c\2c^ 



with V > ^ — 1. 



For a — 1, Gi{t) becomes an exponential random process. The exponential clock permits us to 
derive an interesting interpretation of the probability ( |4. 14| i. Indeed, since 



K. 1 (x) = \ —e 



we get that 



P{X,n{Gi{t)) edx\N{t)=n} = ^[ dww-Hl-wr"^'-^—^e~^ (4.17) 

for n > 0, being the condition > — ^ always satisfied. It means that the conditioned probability of 
Xm{Gi{t)) is equivalent to the distribution of a Laplace random variable with parameter ^ where 

Remark 4.6. The process Xq(Gi (t)) and the occupation time T{t) are linked by the following relation- 
ships 



P{XoiGi{t)) edx\Nit) = 0} 



- f P{T{t) eds} — e-^ds 
TT Jq 2cs 



P{Xo{Gi{t)) e dx\N{t) = 2} ^P{Xi{Gi{t)) e dx\N{t) = 0} 

dx2 



f 1 . 

/ P{T(t) e ds\B{t) < 0} — e^^ds 
Jo 2cs 



as it is possible to derive by ( |4.16[ ) and ( |4.17| l. 



Now, by means of the same steps used in the proof of Theorem 4.2 it in not hard to prove the result 
contained in the next Theorem. Therefore, we omitted the proof. 



Theorem 4.6. For a = 1, we have that 

l-ct ,2 



P{X„(Gi(t)) e dx} = dx r ^-e''^q^{y,t)dy (4.18) 

Jo 2y 

where q"'{y,t) =2p^{y,t), m = 0, 1. 



Theorem 4.6 claims that X,„(t) with exponential time, is equivalent, in distribution, to a Laplace 

^2 

random variable with parameter ^ ^^-^^ . 
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As for the Bessel random time, also in this case we consider a random time obtained iterating I times 
a Gamma random variable. Hence, we indicate the iterated Gamma process with 

tyKi) = G„,(G„,(...(G„,^,(t))...)) 

having density law given by 



CO ,ai foo >OL2 poo fQz+i 



Jo r(ai) Jo T{a2) Jo T{ai+i) 

where aj+i,j = 0, 1, n, is strictly positive. Therefore, by using the same approach contained in the 
proof of Theorem |4.5| the following result yields. 

Theorem 4.7. Given N{t) — n, such that Vm > ^ — I, we have that 



P{Xrn{gi{t))edx\N{t) = n} (4.19) 

r°° +"2 poo y-Qi + i 



dw 



with X € 



Furthermore, at this point also the probability P {Xi,{Qi{t)) e dx\N{t) — 0} immediately follows. 
By setting ai = 1, it is not an hard task to provide the following unconditonal distribution 



P {X,n{Gi{Gc.,{- ■■ {Gc.,^,{t)) ■■■))) e dx} (4.20) 



= dx ^tr-'e-'^'^dh... -tT^^'^e-'^'dti e-^q"^{y,t)dy. 

Jo r(a2) Jo r(a/+i) Jo 2tiy 

Hence, X„i{Gi{Ga2{' • • {Gai^^{t)) •••))) is distributed as a Laplace random variable with parameter 
given by 

e 

G^A---{G^M---)\xUt)y 

In order to conclude the discussion on the random times, we note that and Ga{t) can be mixed 

obtaining a new class of random times. Obviously, R'^{t) and Ga{t) are thought to be mutually indepen- 
dent. Now, we present the following Theorem. 

Theorem 4.8. For n> 1 and Vm > | — 1, we have that 

P {Xra{R''{G^m e dx\N{t) = n} (4.21) 



dx f\ _a V2t fVi\x\\ \^ V2t\x\ 



r(a)S (f,u™ - f + 1) 7o V^c\cV2u}J \^ cVw 

Proof. Since R'^{Ga{t)) has density given by 



r(a)r(^)2i-i Jo z5 
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we can write that 



P {X^{R\G^{t))) e dx\N{t) = n] 

= / P {X^iR^iGo^m e dx\N{t) - n, R\G^{t)) = s] P{R\G^{t)) e ds] 
Jo 

T{v„^ + l)T{v„,) t"dx 



27rr(2w„) r(a)r(|)2 



xj^ ds dz^f^J (c2s2_a;2)-'"-h{|,|<,,}-^e-t^z"-le-*^ 



Therefore, by means of the Corollary 2.1 we are able to explicit the Fourier transform of the previous 
probability distribution as follows 



+00 

e'"^P{X„,{R''{G^{t))) e dx\N{t) = n} 

-00 



r(w™ + !)( — ) ^" , / z°'-^-^e-''-dz I (acs)s'^-''"^-^e-^ds 

t-nv,n+i) r,.-f-i,-..,,f tH f^Y' r,-+^-ie-i^d. 

r(a)r(|)2f-iio ^"^^fc!r(A + t;„ + l)l2J io 

2 \ 



r(a)r(|) Jo 



Jo t^n fc!r(/c + + 1) V 2 ; io 



2 ^ fc=0 

r z'^-'e-'^dzY^ 
r(a)r(f) io ^£^^fc!r(A: + + 



k 



^"^(-+1) ^^"-e-^.E^?^^^^+^™-^l^ ^^^^^^ 



r(a)r(i)r(t,™-f + i)7o t'o ^' V 2' " 2 A 2 

r(a)B(f,,z;„-| + l) io kl \ 2 J J^ 



r(a)B(f,,z;„- 1 + 1) 7o 
By inverting the above characteristic function, we obtain that 

P {X^iR'^iGa^it))) e dx\Nit) = n} 



k=0 

z^'-^e-'^dz [ w^''^{l-wy'--ier''^^dw 



r(a)B(f,,i;„, - 1 + 1) 7o 7o V&^c 

(u = tz) 



-r^ , ^ r wi-\l-wr--Uw /""u^-^-i ^"' "ll^ ' du 

r(a)B(f,,w„-| + l) 7o Jo V2^c 

' dww^ (1 - w) " ^ — ^ if _ - ' 



r(a)B(f,,i;„- 1 + 1) 7o Vc\/2^ 
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□ 

Remark 4.7. For a = 1, one has 



dx , d T d \/2i 



P{X„(i?'^(Gi(t))) e dx\N{t) = n}^ ^— / dww-^-\l-wr--^ 

-D (^2 " ^ 2 + /CVW 

and for d = I, after some calculations, the following distribution yields 

P{X„-,{\BiGim) e dx} - / ^g"(y,i)dy. 



e 



2y 

Remark 4.8. Let X and Ybe two independent random variables distributed as a Gaussian with mean 
zero and variance respectively equal to and 1. It is well-known that V — XY has the following 
density law 

f+-^eSeS 1 /\v\ 



These considerations and Theorem 



4.8 permit us to state that for ct = ^ and d — 1, 



p{x^{\B{G.m)edx\N{t) = n} = —^^ I dwwi-\l-wr--i^K, 



dx . d_,,, _d \/2i _ I \/2i\x\ 

holds for each n. In other words, conditionally on the number of Poisson events during the interval 
[0, t], the law of Xm{\B(G i {t))\) is equivalent to the distribution of the product of a standard Gaussian 

and a Normal random variable, indipendent from previous one, with mean zero and variance ^^=-, with 



5 Some remarks on random motions in higher spaces 

LetX2(i) = {Xi{t),X2{t)) andX^{t) = {Xi{t), X2it), Xs^t), X4{t)) be respectively a planar and 
four-dimensional random flight, then the results presented in the Section |4] can be extended to these 
random processes. 

We also use the follwing notations: 0.2 = (01,02), = (ai, 02, 013, 04), X2 = {xi,X2), X4 = 
{xi,X2,X3,X4). Futhermore, let 1 1 • 1 1 be the euclidean norm and < • , • > the scalar product. As proved 
in Orsingher and De Gregorio (2007b) the characteristic function and the conditional probabilities of the 
random flights in the plane are given by 

while for ^^{t), one has 



\^V'?Mct\\a2\ 
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Therefore, for the planar random flights at Bessel random time, namely {Xi{R'^{t)), X2{R'^{t))), by 
following the same steps of the proof of Theorem |4.2| we have that 

P{Xi(i?'*(i)) e dxi,X2{R'^{t)) e dx2\N{t) = n} (5.1) 

1x2 11^ 



S(f,f -f + 1) Jo 27rtoc2 
P{Xi(G„(t)) e dxi,X2(G„(t)) e dx2|iV(i) = n} (5.2) 

dXidX2 , MJ™-* ^ fA\^2\\\^ / ^112^2 1 



r(s±i)p(f - f + 1) 7o v^cV2cy^y ^ V 

which hold if > f ~ 1 again. We observe that from the probability ( |5.1| l for d = 1, we reobtain the re- 
sult (3.3) showed in Beghin and Orsingher (2009). Similar considerations yield for the four-dimensional 
random flights. 

The random flights in higher spaces have directions uniformly distributed on a multidimensional 
hypersphere. It would be interesting to consider a model, for example in the plane, with a different 
density law with respect to the uniform one. Let us consider a planar random flight with density law 
similar to ( |2.1[ ), for example 

ZTT 

In this case, we obtain a random process describing a motion tending to move in a land of the plane with 
high probability. Therefore, we obtain a random motion with drift, which is persistent along a specific 
portion of the plane. Hence, in order to calculate the characteristic function, we need of the following 
integral 



I exp{iz(acos6' + /3sin6')}sin^ 
Jo 



which it is work out as follows 



p2'K 

i exp{iz(Q; cos 6* + /3sin 61)} sin^ 0^6* 
Jo 

= 5I^r / (acose + f3 sine)'' sin^ Ode 

k=0 ■ 

= inf^Cy^'-' /"cos'^^sin— ^^d^. 



Now, the last integral has to be splitted in two parts: J^^ = + J^^ . Hence, by performing a change 
of variable 9' = 6 — t: in the second integral, we observe by ( |3.1[ ) that the previous sum is equal to zero 
if k and r are odd. By splitting the first integral on (0, |) and after a change of variable analogous to the 
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previous one, we can write that 

I ey:p{iz{acos9 + (3 sm0)}siv? 9d6 
Jo 



fc=0 ^ ' r=0 

OO 9h 



(2fc)!^V2r; " r(fc + 2) 

1).^" f p\ ,,_,, 2-^^r(2r)r(2(fc-r+l)) 
it; (2fc)!^oW r(fc + 2)r(r)r(fc-r + l) 

^ r(fc + 2)22^^4^^2rr(r)r(fc-r+l) ^ 

""ho (fc + l)!2^''' \^2(r!)(fc-r)! ^^r!(fc-r-l)!J 



= 27r. 



/ 9 ^9 Y^—^J2{zVa^ + l3^) 



Then, the characteristic function becomes 



'|e*<S2.X2(*)>|iv(t) = „| (5.3) 



,1 I c(sj- - Sj_i)||a2 

From ( |5.3| l emerges that as expected an asimmetry is introduced by f{0), because the particle will tend 
to maintains the same direction. Moreover, the inversion of the characteristic function is quite difficult, 
therefore seem to be not possible to obtain the explicit probability distribution of X2 at time t by means 
of this approach. 
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